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Abstract 



Within a Two-Higgs-Doublet Model (2HDM) a search for a hght Higgs boson in the mass range 
of 4-12 GeV has been performed in the Yukawa process e"'"e~^ bbA/h — >■ bbr+r", using the 
data collected by the OPAL detector at LEP between 1992 and 1995 in e"'"e~ collisions at about 
91 GeV centre-of-mass energy. A likelihood selection is applied to separate background and 
signal. The number of observed events is in good agreement with the expected background. 
Within a CP-conserving 2HDM type II model the cross-section for Yukawa production depends 
on = I tan/3| and = | sin a/ cos/3| for the production of the CP-odd A and the CP-even 
h, respectively, where tan /? is the ratio of the vacuum expectation values of the Higgs doublets 
and a is the mixing angle between the neutral CP-even Higgs bosons. From our data 95% C.L. 
upper limits are derived for within the range of 8.5 to 13.6 and for i\ between 8.2 to 13.7, 
depending on the mass of the Higgs boson, assuming a branching fraction into t^t~ of 100%. 
An interpretation of the limits within a 2HDM type II model with Standard Model particle 
content is given. These results impose constraints on several models that have been proposed 
to explain the recent BNL measurement of the muon anomalous magnetic moment. 
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1 Introduction 



The search for the last missing particle predicted by the Standard Model, the Higgs boson, 
is one of the main topics in high energy physics at LEP. The Standard Model (SM) has been 
verified to a very high degree of precision but no Higgs boson has yet been discovered. Current 
searches at LEP in the Standard Model scenario [|l| exclude Higgs bosons with masses mno 
below 114.1 GeV at the 95% confidence level. Many proposed models extend the SM while 
preserving the good agreement with experimental data. A minimal extension of the SM, the 
Two-Higgs-Doublet Model (2HDM), has this property. In non-supersymmetric 2HDMs, Higgs 
bosons with small masses still cannot be excluded [0, ^. The analysis presented here, a search 
for A and h in the mass range 4-12 GeV, provides new constraints on the parameter space of 
these models. 

In the 2HDM, two complex Higgs doublets are introduced to generate the mass of the 
fermions: 




One constraint for extended Higgs sectors is the experimental observation that the value of 
p = mw±^/(mz^ cos^ ^w) ~ 1- The condition of p ~ 1 is automatically met by models with 
only Higgs doublets. 

There are several possible patterns for how the new fields may couple to fermions. Four 
types of 2HDMs are theoretically considered in order to avoid introducing flavour changing 
neutral currents (FCNC) The four types differ in the way the two Higgs fields $i and $2 
couple to fermions (see Table 0). 



couples to 


type I 


type II 


type HI 


type IV 


down-type leptons 


$2 


$1 


$2 


$1 


up-type quarks 


$2 


$2 


$2 


$2 


down-type quarks 


$2 


$1 


$1 


$2 



Table 1: The couplings of the Higgs fields according to the four types of 2HDMs. 

In this analysis the CP-conserving 2HDM type II is considered where $1 couples to the 
down-type and the $2 couples to the up-type matter fields. This 2HDM predicts five physical 
Higgs bosons: two charged (H^), two CP-even (H°, h) and one CP-odd (A) Higgs boson (the 
neutral Higgs bosons are often referred to as scalar and pseudoscalar, respectively). The Higgs 
sector of the general CP-conserving 2HDM has six free physical parameters [§] and can be 
parametrized by the masses mH±, rnno, mh and mp^ and two dimensionless parameters a and 
tan/3, with a being the mixing angle in the CP-even neutral Higgs sector, and tan/3 = ^2/^1 
being the ratio of the two vacuum expectation values V2 and Vi. 

The Higgs sector in the Minimal Supersymmetric Standard Model (MSSM) is a special case 
of a 2HDM type II model in which, due to relations imposed by supersymmetry, only two free 



4 



parameters remain (eg. tan (3 and a) at tree level in the Higgs sector. Direct searches at LEP 
in the MSSM scenario have set limits on the Higgs mass and its parameter space. For example 
in a parameter scan where the parameter combination was chosen to give a maximal predicted 
mass mh, 95% CL. limits have been set on the masses mh and larger than 88.4 GeV for 
tan/3 > 0.4 §. 



2 Higgs Boson Production in e+e Collisions 

There are three main processes for Higgs production at Born level within the energy range 
covered by LEP, namely the Higgsstrahlung, associated production and Yukawa processes shown 
in Figure |l|, of which only the first process is of importance in the Standard Model. In the 
Standard Model the Yukawa process is suppressed by the factor {mf/m\i^), and the associated 
production of Higgs bosons is nonexistent. The cross-sections of the 2HDM are closely related 




Figure 1: The three Born level Higgs production processes in a 2HDM at LEP: The Hig- 
gsstrahlung process a), associated production b) and the Yukawa process c). 

to similar processes of the Standard Model: 

a(e+e- ^ Z ^ hZ*) = asuie+e' Z HsmZ*) sm\(3 - a) (2) 

cT(e+e- ^ Z ^ hA) = (TsM(e+e- -> Z ^ z/z/) cos=^(/? - a)\l (3) 
where A = (1 — /«h — ka)^ — 4/th^KA^ being a phase space factor, with = mf/mz^ @|. 

In recent years there have been searches for the Standard Model Higgs boson as well as 
for MSSM Higgs bosons by all four LEP experiments The interpretation of the flavour 
independent Higgs search within the 2HDM in the mass range below approximately 40 GeV 
requires sin^(/3 — a) to be less than 0.2 ||^. For a sufficiently small sin^(/5 — a), the h produced 
through the Higgsstrahlung process can not be seen in the data collected at LEP due to the 
cross-section suppression factor (Eqn. 2). Associated production, if kinematically allowed, 
would be the dominant process for Higgs boson production. On the other hand a light Higgs 
boson which is produced only in the Yukawa process (A or h whichever is lighter) could have 
escaped discovery. Under the assumption that the Higgsstrahlung process is suppressed and 
associated production is kinematically forbidden, mA + > the Yukawa process becomes 
the dominant process for Higgs production at LEP. This scenario can easily be realised in the 
general 2HDM since its parameters are not constrained theoretically In this analysis we 
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concentrate on Higgs masses below the 2mb threshold since previous analyses [0 are insensitive 
to such light Higgs bosons. We will also constrain ourselves to Higgs masses above the 2mr 
threshold. 



The cross-section of the Yukawa process 



CTYukawa OC TTl^ 



(4) 



is proportional to the squared fermion mass, , the colour factor, Nc, of the emitting fermion, 
and an enhancement factor , which describes the coupling between the Higgs boson and the 
emitting fermion (see Table H) . 



Higgs Type 


Down Type Fermions 


Up Type Fermions 


A 
h 


= —sina/ cos P 


e = l/tan/? 
= cos a/ sin/5 



Table 2: The enhancement factor depending on the type of the Higgs boson and the emitting 
fermion. 



The range of for which a detectable signal would be produced can be divided into two 
parts: 

1. < 1 (implies C,^ > 1): An up-type quark pair cc radiates a Higgs boson decaying into 
cc. 

2. > 1: A down-type quark pair bb radiates a Higgs boson decaying into r+r~. 

In this analysis, only the bbr"'"r~ final state is considered, since background suppression 
can be performed more efficiently using the clear signature of bottom decay, together with the 
missing energy and low-multiplicity signature of r decays. In Figure ^ the respective numbers 
of expected events are shown as a function of for the example of a CP-even Higgs with a 
mass of 4 GeV. 



3 Search in the bbr+r Channel 
3.1 Data and Monte Carlo Samples 



The present analysis is based on data collected with the OPAL detector |]T0[ during the years 
1992-1995, taken at centre-of-mass energies close to the Z peak. Although the integrated 
luminosity at LEP recorded at ^/s between 130 and 208 GeV is about a factor five higher 
than the luminosity around the Z peak, the number of produced b quarks is about a factor of 
100 smaller in the higher-energy data. Therefore we only use data recorded at the Z peak at 
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centre-of-mass energies near 91 GeV to search for Yukawa production of Higgs bosons. The 
data collected at off-peak energies in the range between 89 GeV and 93 GeV are included to 
increase the available statistics. Since the characteristics of signal events do not depend on 
the precise center-of-mass energy, and since the background is completely dominated by SM 
hadronic Z decays, we can treat the data as if all were taken on-peak by rescaling the off-peak 
luminosity appropriately. For this we count the number of observed hadronic events [jll| in the 
data. This number, efficiency corrected and divided by the peak cross-section for hadronic Z 
decays of 30.45 nb as measured by OPAL ||12[, then yields an effective on-peak luminosity of 



113.1 pb [IT3, from which about 12% is contributed by the data taken off-peak. 



We consider three types of background classes, with the full detector response simulated as 
using the following generators. 



described in [14 



1. Two-photon background, generated with Vermaseren 1.01 |15| and PHOJET [16 



Two-photon production of hadronic final states is characterized by little visible en- 
ergy in the detector. Usually, the e"*" and e~ escape undetected close to the beam axis, 
causing only a small amount of transverse momentum. The relatively large cross- 
section of this process nevertheless makes us consider it as a potential background. 
The Monte Carlo sample of about two million two-photon events corresponds to a 
luminosity of about four times the data luminosity. 
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2. Four-fermion background, generated with FERMISV [IT^ and grc4f [ITB 



The four-fermion background at LEP 1 mainly originates from Initial State Pair 
(ISP) and Final State Pair (FSP) radiation diagrams. This background class can 
be divided into two subclasses, the first containing four lepton final states and the 
second having qqff final states. The first subclass is eliminated by requiring the 
event to pass a general multihadronic selection [0. The second class contains a 
small amount of irreducible e"'"e^^ hhi^i^ background with a charged lepton pair, 
mainly from FSP radiation. The Monte Carlo sample of about 25000 four-fermion 
events corresponds to a luminosity of about eight times the data luminosity. 



3. qq background, generated with Jetset 7.4 



This class of background consists of events of the type e"'"e~— > Z/7 — > qq. Events 
with gluons radiated off the quarks, especially e"'"e~— > Z/7 ^ bbg(g), are very likely 
to have signal characteristics and represent the main background in this analysis. 
The Monte Carlo sample of about seven million processed qq events corresponds to 
a luminosity about two times the data luminosity. Generated Monte Carlo qq back- 
ground samples for different detector setups are weighted according to the respective 
luminosity for the data. 



Eighteen signal samples of 10000 events each with masses of mh/A=4-12 GeV in one GeV 
steps were generated using a newly written Monte Carlo program based on ^j. The hadroni- 



sation is done with JETSET version 7.408 [19 together with OPAL specific modifications [2C 
The decay of the tau leptons is simulated with the tau decay library TAUOLA [^. These 



signal events are subjected to the same reconstruction and event selection as the real data. 



3.2 Analysis Tools 

In calculating the visible energies and momenta. Ems and Pvis, of individual jets and of the 
total event, "energy-flow objects" are formed from the charged tracks and calorimeter clusters 



22 1 . To avoid double counting, the energies expected to be deposited in the calorimeters by 



the charged particles are subtracted from the energies of the associated calorimeter clusters. 

In order to identify jets containing b hadron decays, three independent techniques using 
lifetime, high-p^ lepton characteristics and kinematic information are used. Artificial Neural 
Networks (ANN's) have been trained to combine several lifetime-sensitive tagging variables and 
kinematic variables. For each jet, the outputs of the lifetime ANN, the kinematic ANN, and 
the lepton tag are combined into a likelihood variable B which discriminates b-fiavoured jets 
from c-fiavoured jets and light quark jets M. 
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3.3 Properties of Yukawa Production 



One of the properties of Yukawa production is the hard energy spectrum of the emitted Higgs 
bosons. In Figure ^ the energy distributions of the A and h are shown. This leads to a high 
boost of the decay products of the Higgs particles and results in a small angle between the 
two r's in the detector. Consequently, the tracks of the decay products of the two r's can 
be reconstructed in one single low-multiplicity jet. In addition, two high-multiplicity jets are 
expected to be associated with the b quarks. Thus the expected topology of a signal event 
has a three jet signature. A simulated signal event is shown in Figure ^ In this analysis each 
event is forced into three jets using the Durham algorithm [Q. The jets are sorted according to 
their multiplicity, assuming that the jet with the lowest charged particle multiplicity (denoted 
'Jet (3)') contains the two r leptons from the Higgs decay. 




Signal m^=4 GeV 
Signal m^=4 GeV 



I .... I .... I .... I .... I .... I .... I .... I 



10 15 20 25 30 35 40 45 50 

Higgs Energy [GeV] 



Figure 3: The energy distribution of the Higgs boson with mh/A = 4 GeV from 10000 Monte 
Carlo signal events at generator level. 



3.4 Event Selection 



The event selection consists of two parts, a preselection and a subsequent likelihood selection. 
Since the unknown Higgs mass enters in the properties of the likelihood variables, the selection 
was performed and optimized separately for each of the nine simulated signal mass hypotheses. 

The preselection consists of the eight cuts described below (see Figure 
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Figure 4: A characteristic simulated signal event e^e~^ bbA/h — > hhr^r^ reconstructed in 
the OPAL detector. Jet(l) and Jet(2) contain the tracks of the hadronizedh quarks and Jet(3) 
contains the decay products of the Higgs boson. The missing momentum vector (dark grey 
arrow) points along the Jet(3) axis due to the undetected neutrinos in the decay of the t's. 



0. General hadronic event selection as described in IT 



1. Eyis < 90 GeV. Due to the neutrinos from r decays the signal events have missing energy 
in the detector. This cut on the visible energy primarily suppresses qq background. 

2. Pt{miss) > 3 GeV. The cut on the transverse missing momentum in the event is addi- 
tionally introduced to reduce two-photon background. 

3. To suppress two-photon events further, we require the missing momentum vector to have 
I cos (beam axis, Pmiss)| < 0.95. 



4. A two dimensional cut 3 ■ C + log(?/32) > —4.5 on the event shape variables, C-value P3 



and ys2 IHI, is introduced to suppress two-jet like events. 

5. A large fraction of the missing momentum in the event is due to the undetectable neutrinos 
of the decaying r's in Jet (3). Therefore a cut on the angle between Jet (3) and Pmiss is 
introduced: cos(Jet(3), Pmiss) > —0.6. 

6. We require at least one identified electron or muon in Jet(3). Electrons are identified 



according to [g^ and muons according to [^. This cut is made to reduce qq background 
and to define efficient likelihood variables. The probability of two r decays containing at 
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Cuts 


Data 


qq 


four-fermion 


two-photon 


(0) 


338 xlO^ 


338 xlO^ 


887 


346 


(1) 


196 xlO^ 


198 xlO^ 


545 


236 


(2) 


131 xlO^ 


128 xlO^ 


357 


123 


(3) 


125 xlO^ 


123 xlO^ 


330 


54 


(4) 


622 xlO^ 


599 xlO^ 


198 


43 


(5) 


400 xlO^ 


383 xlO^ 


125 


32 


(6) 


292 xlO^ 


307 xlO^ 


39 


5 


(7) 


142 xlO^ 


141 xlO^ 


22 


1 



Table 3: The number of events selected and expected in the preselection. The three categories 
of background have been normalised to the data luminosity. 



least one charged lepton is approximately 60%. To ensure the correct efficiency modelling 
in the Monte Carlo a further cut on the lepton momentum, Pi > 2 GeV, is added. 

7. 2 < N track < 4 in Jet (3). This cut on the charged multiplicity of Jet (3) is introduced in 
order to optimize the ratio of signal over square root of background after the likelihood 
selection. 

As shown in Table the observed number of data events and the expected background 
agree well at each step of the preselection. Of the background passing the general multihadronic 
selection (cut 0) about 2% of the four-fermion events, 0.4% of the qq events and 0.3% of the 
two-photon events remain after the complete preselection, with the qq background clearly being 
the dominant contribution. The signal selection efficiency, only weakly dependent of the Higgs 
boson mass, ranges from 11% to 17%, depending on the signal sample, as detailed for some 
typical masses in Table ^ 



Cuts 


Efficiency (%) 


Efficiency (%) 


Efficiency (%) 


Efficiency (%) 




mA=4 GeV 


mA=10 GeV 


mh=4 GeV 


mh=10 GeV 


(1) 


92.0 


92.7 


90.2 


91.6 


(2) 


80.6 


83.0 


76.7 


81.4 


(3) 


78.7 


80.8 


74.9 


79.4 


(4) 


64.5 


73.9 


54.9 


71.2 


(5) 


54.8 


65.4 


43.8 


60.7 


(6) 


19.2 


22.2 


14.2 


19.3 


(7) 


16.2 


17.4 


11.2 


14.7 



Table 4: The efficiency of the preselection for selected Higgs masses. 



From reference histograms of eight variables, listed below, we define signal likelihood selec- 
tions for each mass hypothesis, both for the h and the A Higgs boson. Due to the overwhelming 
dominance of the qq background, we use a single inclusive background class. The eight reference 
histograms are (see Figure |^) : 
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1. Bi. The 'b-ness' of tlie jet witli tlie liigliest multiplicity. This value is defined as 



Lb + L, + L 



uds 



(5) 



Here Lf,, Lc and Luds are likelihood values for bottom, charm and light flavour jets re- 
spectively [0]. 

2. i?2. The 'b-ness' value of the jet with the second highest multiplicity. 

3. Mms- The measured invariant mass of the event. 

4. The C- value of the event ||2^ . 



5. (Pi -|- P2)/ Ejet in Jet(3). The sum of the momenta of the two tracks in Jet(3) with the 
highest momentum divided by the measured energy in Jet (3). 

6. Pt- The transverse momentum of the event, with respect to the beam axis. 

7. log(?/32). The logarithm of the 1/32 value of the event p4 |. 



8. cos(Jet(l), Jet(2)). Cosine of the angle between the jet with the highest multiplicity 
(Jet(l)) and the one with the second highest multiplicity (Jet(2)). 



The likelihood distributions are shown in Figure |^. The likelihood cuts are determined sep- 
arately for each Higgs type and mass hypothesis in a compromise to achieve a good expected 
limit, calculated with Monte Carlo experiments, and a smooth behavior of efficiency and ex- 
pected number of backgrounds as a function of the Higgs mass. After the likelihood cut the 
data are in good agreement with background Monte Carlo simulation (see Table ^and Table |^). 



Mass(A) 
(GeV) 


LH cut 


Data 
Events 


Total background 
Events 


Efficiency 

% 


4 


0.985 


12 


14.9±2.9±1.6 


3.4±0.2±0.1 


5 


0.985 


17 


16.0±3.1±1.7 


3.5±0.2±0.1 


6 


0.985 


13 


16.6±3.1±1.7 


3.7±0.2±0.1 


7 


0.987 


14 


13.8±2.8±1.4 


3.8±0.2±0.2 


8 


0.990 


13 


17.8±3.3±1.9 


3.7±0.2±0.1 


9 


0.990 


11 


15.1±3.0±1.6 


3.9±0.2±0.2 


10 


0.990 


11 


17.6±3.2±1.8 


3.8±0.2±0.2 


11 


0.990 


13 


19.1±3.3±2.0 


4.0±0.2±0.2 


12 


0.992 


13 


18.3±3.3±1.9 


3.9±0.2±0.2 



Table 5: The number of selected A candidate events after the likelihood cut. 
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Mass(h) 
(GeV) 


LH cut 


Data 
Events 


Total background 
Events 


Efficiency 

% 


4 


0.965 


38 


46.9±5.1±4.8 


2.8±0.2±0.1 


5 


0.970 


41 


36.9±4.6±3.8 


3.1±0.2±0.1 


6 


0.975 


26 


29.0±4.1±3.0 


2.8±0.2±0.1 


7 


0.980 


22 


20.5±3.4±2.1 


3.2±0.2±0.1 


8 


0.985 


18 


12.7±2.7±1.3 


3.1±0.2±0.1 


9 


0.987 


16 


15.3±3.0±1.6 


3.2±0.2±0.1 


10 


0.985 


15 


17.9±3.2±1.8 


3.0±0.2±0.1 


11 


0.990 


11 


10.1±2.4±1.0 


3.3±0.2±0.1 


12 


0.993 


9 


9.9±2.4±1.0 


3.3±0.2±0.1 



Table 6: The number of selected h candidate events after the likelihood cut. 



4 Systematic Studies 



In order to estimate systematic uncertainties, the following sources of uncertainty are varied 
in the Monte Carlo, one at a time, and the analysis is redone with the adjusted samples to 
calculate the difference in the background expectation and the efficiency with respect to the 
standard analysis. For variations done in two directions, with respect to the central value, 
the largest deviation in expected background and efficiency after re-analysis is taken as the 
systematic uncertainty. The following sources of uncertainty were studied: 



The systematic error contributed by the b-tagging due to track resolution modelling. 
This was estimated with a variation of the track parameters 0o i 5%, do ± 5%, zq ± 
10% as described in [^. This variation particularly influences the values of track- and 



secondary vertex significances {d/ad and l/cri) which are important variables for the b- 
tagging procedure. 



The uncertainty in the b-quark fragmentation function |2^ is estimated by varying eb by 
±25% around a central value of eb = 0.0038 A smaller (higher) value of eb corresponds 
to a harder(softer) B-hadron spectrum. 

Description of the kinematic variables used in the likelihood selection. The kinematic 
likelihood variables in the Monte Carlo were shifted to match the mean of the data. After 
the shift on the variable the selection is reapplied and the deviations in the expected 
background and efficiency are taken as systematic errors. 



The B-hadron charged decay multiplicity uncertainty |30 



The systematic deviations are calculated separately for the two different hypotheses on 
Higgs type (A and h) after a likelihood cut at 0.8 to ensure that statistical contributions to the 
estimated systematic errors are minimal. 
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A Selection 


li Selection 




Signal eff. 


qq background 


Signal eff. 


qq background 


Variation of 


% 


Events 


% 


Events 


Track parameters 


0.2 


4.1 


2.4 


1.7 




2.7 


8.3 


1.7 


9.2 


B multiplicity 


1.8 


1.8 


1.8 


1.8 




0.3 


2.6 


0.9 


1.8 


(Pi + P2)/E 


0.3 


0.3 


0.1 


0.2 


C-value 


0.6 


3.4 


0.8 


2.7 


Pt 


0.1 


0.6 


0.1 


0.5 


log(y32) 


0.5 


0.8 


0.2 


0.2 


cos(Jet(l), Jet(2)) 


1.7 


1.1 


0.3 


1.1 


r correlations 


1.3 


1.5 


1.3 


1.5 


Total 


4.0 


10.5 


3.9 


10.3 



Table 7: Relative change of the signal efficiency and number of expected background events for 
the selection of A and h at a likelihood cut, L > 0.8. The total value expresses the quadratically 
added contributions. 



Furthermore the analysis for one typical Higgs hypothesis (mh=5 GeV) was redone under 
the assumption of 100% correlated helicity states of the r's. The deviation of 1.3% for signal 
and 1.5% for background expectation were included as a source of systematic uncertainty. The 
contributions to the systematic uncertainty are broken down for the analyses in Table |^. 

Adding in quadrature the statistical uncertainties and the uncertainties from the above 
sources yields the total errors, listed in Table ^ and Table ^, on the selection efficiency and 
background rates for all nine masses of A and h. The uncertainty due to the b-fragmentation 
function is the highest contribution. 



5 Limits on ^^/^ in the 2HDM Type II 

The results of the selection with systematic and statistical errors are listed in Table |^ for the 
A and in Table ^ for the h. 

2HDM Type II limits are determined for the cross-section of the process e^e^— Z— bbr^r^. 
Due to the dependence of the cross-section (Eqn. 4) on the enhancement factor C,d^^, a limit on 
^ is obtained by calculating A^95, the 95% C.L. upper limit on the rate of accepted signal events 



in the data, according to and adjusting for the efficiency, cross section, and luminosity. 



The total error on the efficiency and on the background estimation are convoluted into the 



limit according to [32|. For the mass points between the measured ones, in one GeV steps from 



4 to 12 GeV, the limit is linearly interpolated from the two neighbouring measurements. Com- 
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paring limits obtained with analyses optimized to neighbouring mass points, the uncertainty of 
the interpolation was estimated to be less than 0.5 units in ^. 



Type 


^h/A GeV 


4 


5 


6 


7 


8 


9 


10 


11 


12 


A 


expected 


9.5 


10.3 


11.0 


11.1 


12.9 


13.4 


14.4 


15.5 


16.6 


A 


observed 


8.5 


11.0 


9.6 


11.5 


10.7 


11.0 


11.3 


12.3 


13.6 


h 


expected 


10.2 


10.3 


11.3 


10.8 


10.5 


12.1 


13.7 


11.9 


12.9 


h 


observed 


8.2 


11.8 


10.4 


11.8 


13.7 


12.6 


12.6 


12.7 


12.9 



Table 8: The upper limit on ^ in a 2HDM Type II model for masses muiggs = 4 GeV to 
f^Higgs = 12 GeV at 95% C.L. calculated assuming 100% branching ratio of the A (upper part) 
and h (lower part) into t^t^ . 



Assuming a branching ratio of 100% for Higgs boson decays into r+r the limits on are 
shown in Figure ^(a) for A production and in Figure ^(b) for h production, and are summarised 
in Table g 

In a 2HDM model with Standard Model particle content, the Higgs branching ratio into 
T^T~ for 10 is about 85% for Higgs masses between 4 and 9.4 GeV. In the mass range 
from 9.4 Gev to 11.0 GeV the branching ratios are very much influenced by mixing of the Higgs 
bosons h and A with bb bound states with the same quantum numbers (see Table |^). We have 
therefore calculated the branching ratios of the Higgs into r+r^ according to reference . 
The limits derived within this model are shown in Figure |^(a) and (b) for CP-odd and CP-even 
Higgs production, respectively. The observed structure, in particular at higher Higgs boson 
masses, is a consequence of the behaviour of the branching ratio in this particular model. 



state 


Tj Mass 


Xo Mass 


n 


GeV 


GeV 


1 


9.412 


9.860 


2 


9.992 


10.235 


3 


10.340 




4 


10.570 




5 


10.846 




6 


11.014 





Table 9: The mass of the rj states assumed to mix with the pseudoscalar A and of the xo states 
assumed to mix with the scalar h, taken from j^/ . 
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6 Implications for the Muon Anomalous Magnetic Mo- 
ment 



Recent measurements of the anomalous magnetic moment, ct^ = ^{g — 2)^, of the muon have 
given a result which deviates from the Standard Model expectation by 
sponding to about 2.6 standard deviations [B^ 
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400 X 10' corre- 
Depending on the estimation of the hadronic 
contribution to the muon anomalous magnetic moment, the 90% CL. ranges for the contribu- 
tion of New Physics 5a^(NP) are: 



215 X 10"^' < 5a^(NP) < 637 x 10"^' [35] 
170 X 10~^^ < 5a^(NP) < 690 x 10^ 
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112 X 10"^' < 5a^(NP) < 573 x 10"^' [p7 



(6) 
(7) 
(8) 



Light Higgs bosons A and h could form a part of via loop diagrams. A one-loop calcula- 
tion [0 predicts positive contributions 5a^(h) > for the h, and negative contributions 5a^(A) 
< for the A. The two-loop terms, due to the stronger coupling of the Higgs fields to loops 
with heavy quarks, turn out to be larger in magnitude than the one-loop terms, and of opposite 
sign , giving a total positive contribution 6aj/{A) > for the A, as shown with indicated 



isolines in Figure ^ (a) [^. However, the two- loop terms gives a total negative contribution 



Saj/{h) < for the h, thus suggesting that the h can not account for the BNL observation. 
We show in Figure ^ (b) only the isolines of the contribution from the earlier one-loop calcula- 
tion which resulted in a positive value of (5a^(h). Our data exclude positive contributions 
(5a^(h) > 100 X 10^^^ for h masses between 4.0 and 10.7 GeV at the one-loop level, and 6a^J{A) 
> 100 X 10^^^ for A masses betwen 4.0 and 9.9 GeV at the two-loop level. Similar limits have 



been derived from radiative T decays [^, Q for h/A masses lighter than about 8 GeV, with, 
however, large QCD uncertainies. 



In reference and |^ the authors have suggested that a light Higgs boson could fully 
account for the observed deviation of the measured {g — 2)^ from the Standard Model expec- 
tation. In a scenario without contributions of other new particles, eg. gauginos, and assuming 
that either h or A is heavy enough to render associated Ah production inaccessible at LEP, 
only the lighter of the two Higgs bosons would sizeably contribute. For a light h, one has to 
assume in addition sin(/5 — a) ~ to explain its non-observation in the Standard Model search 
for the Higgsstrahlung process. The experimental results of this analysis can be interpreted in 
such a scenario and would rule out (using the one-loop calculation [^) a light h in the mass 



range of 4-10.7 GeV, and (using the two-loop calculation a light A in the mass range from 
4-9.9 GeV as the only source of the discrepancy in the {g — 2)^ measurement for all three 90% 
C.L. ranges listed above. 
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7 Conclusion 



The Yukawa production of light neutral Higgs bosons in the channel e"'"e~— > bbA/h— >^ bbr"'"r~ is 
studied. The search presented here, based on data collected by OPAL at ^/s ~ mz in the years 
1992 to 1995, has not revealed any significant excess over the expected background. Limits on 
the Yukawa production of a light Higgs with masses in the range of 4 GeV to 12 GeV have 
been set at 95% C.L. New limits on the parameters = |tan/?| and = |sina/ cos/3| are 
presented for A and h production, respectively. Assuming a branching ratio to r+r^ of 100% 
, upper limits for can be set within the range of 8.5 to 13.6 and for between 8.2 to 13.7, 
depending on the mass of the Higgs boson. In a 2HDM type II model with Standard Model 
particle content similar limits are obtained up to masses of 9.4 GeV. Above 9.4 GeV the mixing 
of the Higgs bosons h and A with bb bound states with the same quantum numbers, especially 
for the CP-odd A, results in weaker limits in certain mass ranges. The experimental result of 
this analysis restricts the contribution to the anomalous magnetic moment of the muon for a 
light h in the mass range of 4-10.7 GeV (using one-loop calculation I^H]) and for a light A in the 
mass range from 4-9.9 GeV (using the two-loop calculation of [^) to (5a^(Higgs) < 100 x 10~^^ 
at the 95% C.L. 
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Figure 5: The preselection variables and their impact on processed data. The points with error 
bars are data, while the solid histogram is the simulation of the qq background normalised to 
the recorded luminosity. The dashed line represents a simulated signal of a scalar Higgs with 
rrih = 4 GeV scaled arbitrarily for better visibility. The arrows indicate the cuts made on the 
variables. 
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Figure 6: Distributions used in the likelihood selection for preselected events and two hypothetical 
h Higgs masses. The points with error bars are data. The solid line is simulated background 
normalised to the recorded data. The dashed (dotted) line is a simulated Higgs boson h at a 
mass mHiggs = 4 GeV (mmggs = 10 GeV^ scaled arbitrarily for better visibility. 
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Figure 7: Likelihood distributions for the selection of a CP-odd Higgs A with masses of 4 and 
10 GeV and for a CP-even Higgs h with masses of 4 o,nd 10 GeV. The points with error bars 
are data. The solid line is simulated background normalised to the recorded luminosity. The 
dashed line is a simulated Higgs boson scaled arbitrarily for better visibility. 
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Figure 8: Excluded values of o^t 95% C.L. (dark grey region) in the 2HDM type II model for 
the Yukawa production of a CP-odd A (upper plot) and for the CP-evenh (lower plot) assuming 
the branching fraction of the Higgs boson into r^r^to be 100%. The expected (diamonds) and 
observed (stars) limits have been calculated at specific mass points (4^12 GeV in one GeV steps) 
and linearly interpolated in between. 
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OPAL 




Figure 9: Excluded values of at 95% C.L. in the 2HDM type II model with Standard Model 
particle content for the Yukawa production of a CP-odd A (upper plot) and for the CP-even 
h (lower plot) with the mixing to bb hound states taken into account. The structure results 
entirely from the theoretically suggested branching ratios /|^/. The expected (diamonds) and 
observed (stars) limits have been calculated at specific mass points (4^12 GeV in one GeV 
steps) and linearly interpolated in between. The dotted lines are the contours of the predicted 
Higgs contribution (one-loop [^ / for the CP-even and two-loop / p^ , |^ for the CP-odd 
Higgs) to the muon anomalous magnetic moment, 5a^(Higgs) (in units oflQ^^"^). 
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